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ABSTRACT 

We present the recently updated and expanded MASSCLEANco/or^, a database of 70 million Monte 
Carlo models selected to match the properties (metallicity, ages and masses) of stellar clusters found in the 
Large Magellanic Cloud (LMC). This database shows the rather extreme and non-Guassian distribution of 
integrated colors and magnitudes expected with different cluster age and mass and the enormous age de- 
generacy of integrated colors when mass is unknown. This degeneracy could lead to catastrophic failures 
in estimating age with standard SSP models, particularly if most of the clusters are of intermediate or low 
mass, like in the LMC. Utilizing the MASSCLEANc6>/6>r^ database, we have developed MASSCLEANag^, 
a statistical inference package which assigns the most likely age and mass (solved simultaneously) to a 
cluster based only on its integrated broad-band photometric properties. Finally, we use MASSCLEANag^ to 
derive the age and mass of LMC clusters based on integrated photometry alone. First we compare our clus- 
ter ages against those obtained for the same seven clusters using more accurate integrated spectroscopy. 
We find improved agreement with the integrated spectroscopy ages over the original photometric ages. A 
close examination of our results demonstrate the necessity of solving simultaneously for mass and age 
to reduce degeneracies in the cluster ages derived via integrated colors. We then selected an additional 
subset of 30 photometric clusters with previously well constrained ages and independently derive their age 
using the MASSCLEANag^ with the same photometry with very good agreement. The MASSCLEANag^ 
program is freely available under GNU General Public License. 
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1. Introduction 

Stellar clusters are among the most universally- 
useful entities studied by modern day astronomers 
Having a set of singly-aged, same-distance stars 
in an identifiable cluster provides stellar and 
galactic astronomers with powerful diagnostic 
tools to interpret our universe. Of particular in- 
terest is the ability to use the ages and masses 
of constituent stellar clusters to track the star for- 
mation history and chemical evolution of a galaxy, 
the mass function of stellar clusters and their life- 
time vjnthiji^^ 



2QQ5I : iPiatti et aLll^QOl : lLarsenll2QQ9l : iMora et aL 



20091 : IChandar et ail 1201 Oah . However, such 
measures require deriving the ages and masses 
of stellar clusters when the constituent stars 
within those clusters are not individually resolved. 
With distant stellar clusters, one has only the in- 
tegrated properties to work from, typically in ei- 
ther broad band photometry or inte grated spec- 
troscopy, to estimate age s (e.g. ISearle et aL 
Il980l : Ivan den Berahlll98lh . In such cases, we 
rely on sophisticated models of these star sys- 
tems to derive the cluster properties (age, mass, 
chemistry) based on their observed, bulk proper- 
ties. For this reason there exists a long history 
of effort interpreting integrated cluster proper- 
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ties throu gh the use of these sophisticated mod- 
els (e .g. 'Girardi et al.' '1995'; 'Bruzual & Charlo^ 



■1 



2. The MASSCLEANco/ors Database - Now 
Over 70 Million Monte Carlo Simulations 



2003 : LCer vino & Valls-Gabaud 2009; 



1989l:lChios i 1989; Bruzual 2002; Bruzual 201 
Gerviho & Luridiana 2004; Cervino & Luridiana 



2006; Faaiolini et al. 2007; Lancon & Mouhcine 



2009; 



1999 



Fouesneau & Lancon 2010"; 'Brocato et al.J 



We have computed the integrated colors and 
magnitudes as a function of mass and age for 

a simple stellar cluster. Our results - the newest 

2000l:lLancQn & Fouesneau 2009l:lFQuesneau & Lanconl version of the MASSCLEANc^/^r^ database l|Popescu & HansonI 

1201 Oal ) - are based on over 70 mi llion Monte 
Carlo simulations usi ng Padova 2008 JMar go et aL 



Brocato et al. 2000; Cantiello et al. 2003; 
Raimondo et al. 2005; Raimondo 2009; Santos et liTl 
2006; Gonzalez et al. 2004; Gonzalez-Looezlira eTaD 
2005 ; Gonzalez-Lopezlira e t al. 2010; Maiz Apella niz 
20091: iPessev et all l2009l : iPopescu & HansonI 
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Quite recently, we presented our own inde- 
pendent effort to provide models for interpreting 
integrated stellar cluster photometry. Our model, 
MASSCLEA N (MASSive CLuster Evo lution and 
ANalvsiJI - iPopescu & HansonI I2009I ). is rather 
unique in its design compared to other stellar 
cluster models and provides additional insight 
to interpreting presently available data on stel- 
lar clusters. In this paper we first present in 
§2 our recently expanded data-set of 70 mil- 
l ion Monte Carlo niodels, MASSCLEANc(9/(9r^ 
IjPoDescu & Hanson' '201 OaV This database of 
models shows the stochastic and non-Gaussian 
behavior of integrated stellar cluster colors and 
magnitudes as a function of cluster mass and 
age. This, our first, newly completed, extended 
database of models, was selected to match the 
mass, age and metallicity range expected in the 
stellar clusters of the Large Magellanic Cloud 
(LMC). In §3, we demonstrate the utility of this 
extensive database. When used with a statis- 
tical inference code we have designed within 
MASSCLEAN called, MASSCLEANag^, we use 
the database to work backwards from the inte- 
grated magnitudes and colors of real LMC clus- 
ters, to derive their stellar ages and masses. We 
also discuss our new age results as compared to 
previous age determinations based on the same 
photometric data, to provide initial tests on the 
cluster results obtained using MASSCLEANag^. 
Final discussion of our results and conclusions 
are found in §4. 



^ http : //www . physics . uc . edu/~popescu/massclean/ 

MASSCLEAN package is freely available unde GNU General 
Public License 



20081 ) and Geneva (Lejeune & SchaereH I2001I) 
isochrones, and for two metalicities, Z = 0.019 
(solar) and Z = .008 (LM C). The simula tion were 
done using Kroupa IMF l|Kroupal[2003 ) with 0.1 
Mq and 120 Mq mass limits. The age range 
of [6.6,9.5] in loq{aqe/yr) was chosen (similar 
to iPopescu & HansonI I2009I : iPopescu & HansonI 
IgOlOa) to accomodate both Padova and Geneva 
models. The mass range is 200-100,000 M©. An 
average of 5000 clusters were simulated for each 
mass and age. 

In this work we will look closely at a subset of 
the MASSCLEANc6>/or^ database. We present in- 
tegrated magnitudes and colors (My, (U - B)o, 
{B - V)o) representing 50 million Monte Carlo 
simulations using Padova 2008 stellar evolution- 
ary models alone and for the single metallicity, 
Z = .008 ([Fe/H]=-0.6), selected to match the bulk 
properties of the LMC over m ost of the past few 
billion years (jPiatti et al.ll2009l) . 



2.1. 



The Distribution of Integrated Colors 
and Magnitudes as a Function of Mass 
and Age 



In IPopescu & Hansonll2010al . we analysed the 
influence of the stochastic fluctuations in the 
mass function of stellar clusters on their inte- 
grated colors as a function of cluster mass. Using 
13 mass intervals ranging from 500-100,000 M©, 
we showed the mean value of the distribution of 
integrated colors varies with mass. We further 
showed that the dispersion of integrated colors, 
measured by la (standard deviation), increases 
as the value of the cluster mass decreases. In 
this new work we will present not just the sigma 
range and mean, but the full distribution of inte- 
grated colors and magnitudes. We are also using 
a much larger number of Monte Carlo simulations 
and cluster mass values in the 200-100,000 M© 
interval. This allows us to look closely at where 



MASSCLEANag^ - Stellar Cluster Ages from Integrated Colors - 



the distribution of real cluster properties lie in the 
critical color-age and color-color diagrams. 

The distribution of integrated My magnitudes 
as a function of mass and age is presented in the 
upper panel of FigureQ] The corresponding mass 
values are color-coded and presented in the bot- 
tom panel. Sixty-five (65) different mass intervals 
have been explored in the MASSCLEANc6>/6>r^ 
database as part of this demonstration. All 65 
mass intervals are plotted on top of each other, 
making it difficult to see the individual distribu- 
tions. However, to first order, higher mass clus- 
ters have higher absolute magnitudes. My. 

In Figures [2] and [3] we present the distribution 
of {B - V)q and (U - B)q integrated colors as a 
function of mass and age. The mass values are 
represented by the same colors as shown in the 
bottom panel of Figure Q] The variation of inte- 
grated colors computed in the infinite mass limit 
(10^ Mq in our simulations) and consistent with 
standard simple stellar populati on (SS P) models 
based on Padova 2008 (Mari ao et al.l [2008). is 
represented by th e solid white line in both figu res. 
As was shown in lPopescu & Hansorill2010al . the 
dispersion in integrated colors grows larger for 
lower mass clusters, both in the blue and red side. 
As the cluster mass increases, the dispersion de- 
creases and the integrated colors approach the 
infinite mass limit shown in white. While it is valu- 
able to see the range of masses over-layed in a 
single plot, it is difficult to see the exact shape of 
the color distributions underlying each mass in- 
terval. 

2.2. The Shape of the Distribution of Inte- 
grated Colors and Magnitudes 

We have selected a few single values of mass 
for independent plots in Figures SHgl The dis- 
tribution of integrated My magnitudes is pre- 
sented in Figure [H for clusters mass of 200, 
1,000, 5,000, 10,000, 25,000, and 50,000 M©, 
respectively. The colors are the same as the 
ones used in Figure Q] What is also given in 
a black line is the mean value of the distribu- 
tion of integrated My mag nitudes for that mass 
l|Popescu & Hansonir2010ah . Several very critical 
effects are immediately apparent. The distribu- 
tion (range) of observed integrated magnitudes 
for low mass clusters gets very large with lower 



mass clusters. This was already presented in 
Popescu & Hanson|[2010al . However, the shape 
of that dispersion was not presented like we have 
done here. Here we see in Fig. [H particularly 
with 200 and 1 ,000 solar mass clusters, the dis- 
tribution becomes bimodal for the lowest mass 
clusters, especially for young ages (log t < 8.0). 
A cluster with no more than 200 solar masses 
can, at times, have an absolute magnitude typical 
of a cluster 25 - 100 times more massive and at 
the same age. Even if one can derive an accu- 
rate age for their cluster, one can not rely on the 
absolute magnitude alone to derive mass unam- 
biguously. 

In Figures [5] and [6] we present the distribution 
of integrated colors {B - V)o and {U - B)o as 
a function of age for the same select values of 
cluster mass, 200, 1,000, 5,000, 10,000, 25,000, 
and 50,000 M©, respectively. The integrated col- 
ors as a function of age, computed in the infinite 
mass limit (lO^M© in our simulations) is repre- 
sented by the solid black line. Here the structure 
in the values observed start out closely aligned 
with the classical SSP codes predictions (infi- 
nite mass limit) for the high mass clusters, but 
quickly fans out to a rather broad distribution be- 
low 1 0,000 Mo . When one gets below 1 ,000 M© , 
the color range is extreme, both in the red and 
blue colors predicted by our models. Moreover, 
the distribution shovys a clear bimodal distribution 
iPopescu & Hanson] |2010a|-|Lanpon & Moujicine 




ICerv 



In 



no & Luridiana 2004; Faaioli ni et al 
Popescu & Hanson 2010a . we dis 



cussed the increased dispersion in colors as the 
cluster mass decreased, though we did not inves- 
tigate masses as low as shown here. However, 
even clusters with masses of 1 ,000 to as much as 
5,000 Mq show an extraordinary range in colors, 
away from the predicted SSP single expectation 
value with age. 

2.3. The Distribution of Integrated Colors in 
the Color-Color Diagram 

All the data from Figures[2]and[3]are again pre- 
sented in Figure [7] as a {U - B)o vs. {B - V)o 
color-color diagram. And, once again, the varia- 
tion of integrated colors computed in the infinite 
mass limit (as predicted by classical SSP codes) 
is represented by the solid white line. The 65 
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Fig. 1 . — The results from 50 million Monte Carlo simulations created with MASSCLEAN. Models were generated with 
65 different masses and evolved following Padova 2008 evolutionary models. Higher mass clusters are plotted over 
preceded lower mass clusters. This figure gives the (partially correct) impression that Mv monotonically increases with 
cluster mass for all ages. 
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Fig. 2. — The results from the same 50 million Monte Carlo simulations as shown in Figure [Hand using the same 
color coding to represent cluster mass. The range of (B - V)o colors, as a function of mass and age is given. Again, 
higher mass clusters are plotted over preceded lower mass clusters, so covers the real distribution of expected colors 
as a function of mass and age. A white line which lies at the color center, represents the color dependence with age 
computed in the infinite mass limit (lO^M© in our simulations). This figure gives the (correct) impression that the range 
of cluster colors observed as a function of age, is greatly increased for lower mass clusters. 
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Fig. 3. — The same as Figure 2, only tliis figure siiows (U - B)o colors as a function of mass and age. 
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Fig. 4. — Beginning in the upper left panel, the black line shows the mean value of the distribution of integrated Mv 
colors for a 50,000 M© cluster. The red points are MASSCLEAN, Monte Carlo derived observed magnitudes for a 
50,000 Mq cluster, as a function of age. Over one million models are used to derive the red points in this figure. The 
distribution seen about the mean value is relatively small. In the next panel to the right, is given the magnitudes, as 
a function of age, for Monte Carlo derived cluster with mass, 25,000 M©. The black line is the mean value of the 
distribution of integrated My colors for a 25,000 M© cluster. In the following 4 subfigures, we provide lower and lower 
mass simulated clusters, 10,000 (upper right), then 5,000 (lower left), 1,000 (lower center) and finally 200 M© (lower 
right), all the while showing the mean value of the distribution of integrated My colors for the corresponding clusters 
mass. For clusters with masses of only a few thousand solar mass, the observed range becomes extremely large, and 
can even mimic clusters of lower or even considerably higher total mass. 
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Fig. 5. — Like with FigureH we iiave selected six mass examples (50,000, 25,000, 10,000, 5,000, 1000, and 200 M©) 
to show the observed distribution of (B - V)o colors as a function of age, compared to the infinite mass limit (in the 
black line, the same for all six plots). For the higher mass clusters, 50,000 and 25,000 M©, the range of colors is small 
and stays reasonably close to the values computed in the infinite mass limit. However, in the remaining panels one 
sees that the observed values from our simulations can demonstrate rather extreme colors, bluer, but more so very red 
colors at times, for increasingly lower mass clusters. 
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Fig. 6. — This is the same as for Figure[5l only here we show the variation with mass and age of (IJ - B)o colors, as 
derived through Monte Carlo simulations using MASSCLEAN. 
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masses given are color-coded similarly to Figure 
[H Provided the mass of the cluster is greater than 
50,000 Mq (in the pink), the predicted properties 
of our Monte Carlo generated clusters stays rea- 
sonably close to the integrated colors computed 
in the infinite mass limit. For clusters with masses 
less than 10,000 M© (the dark green just beyond 
the yellow), the range of possible colors becomes 
quite large, overlapping heavily with many differ- 
ent aged clusters and causing enormous degen- 
eracy problems. 

The classical method for cluster age determi- 
nation using integrated colors is based mainly on 
the {U - B)o vs. {B - V)o color-color diagram. 
The large dispersion presented in FigureHshows 
the limitation of this method when only the infi- 
nite mass limit is used. To better view what is 
happening in the complex Figure [71 we show an- 
other view of this data in Figure [8| Here, we ap- 
ply an entirely new color scheme to represent our 
Monte Carlo simulations. In Figure [8] we plot the 
color-color diagram for 5 ages (log{age/yr)= 7.00, 
7.50, 8.00, 8.50, 9.00) and all values of mass be- 
tween 200-100,000 Mq. Clusters of all masses 
are shown now color coded based on age. Com- 
bining this with what we have learned from Figure 
[71 Figure [8] shows us that even very young low 
mass clusters can mimic the UBV colors of more 
massive old clusters. The extreme tail of colors 
belongs solely to intermediate and low mass clus- 
ters. It is easy to see that very young clusters 
(log{age/yr)= 7.00, 7.50) could be mistaken for 
1 billion years old (log{age/yr)= 9.00) clusters, 
based on the single SSP relationship in a clas- 
sical UBV diagram. 

The {U - B)o versus {B - V)o color-color di- 
agram is a very important tool in age determina- 
tion of stellar clusters. The high level of dispersion 
shown in Figures [7] and [8] demonstrates the limi- 
tations of this method when only the infinite mass 
limit is used. The degeneracy seems to be so ex- 
treme, one might think all is lost. But in fact, the 
purpose of our extensive database was not just 
to show the range and degeneracy of the colors 
with age and mass (and here we have only lim- 
ited ourselves to a single IMF, chemistry and evo- 
lutionary model), but to see if by producing such 
a database, we might provide a reasonable grid 
of properties from which to work back and derive 
cluster ages and masses. 



2.4. Mass and Age Are Degenerate Quanti- 
ties in UBV Colors 

Figures [H through [8] demonstrate the very dif- 
ficult condition one faces in deriving unique ages 
of stellar clusters when mass is not given full con- 
sideration. In creating the MASSCLEANco/(9r^ 
database, we have created a distribution func- 
tion of observed photometric properties with age 
and mass. Our figures demonstrate just the 
UBV bands, however, our MASSCLEANc(9/or^ 
database includes all photometric bands presently 
included in the Padova 2008 model release: UB- 
VRIJHK. While our figures seem to indicate there 
is little to match on, in fact, when simultaneously 
solving all three unique values of UBV (U - B, 
B -V and Mv) , there will be optimized, or most 
likely solutions when one searches in both age 
and mass grid space. If additional photometric 
bands are available for a cluster, this can be fur- 
ther used to constrain the cluster properties. The 
module which will provide such a search we call 
MASSCLEANag^. 

3. Age Determination for Stellar Clusters - 
MASSCLEANa^^ 

The newest addition to t he massclean package 
([Popescu & Hansonl 120091) is MASSCLEANagfl 
a program whic h uses the MASSCLEANc6>/6>r^ 
database ( Popescu & Hanson 2010a) to deter- 
mine the ages of stellar clusters. 

Our program, MASSCLEANag^, is based on a 
method analogous to x^ minimizatiorEl. Based 
on the observational data, integrated magnitude 
and integrated colors, it searches for the mini- 
mum hyper-radius \n the integrated magnitude, in- 
tegrated colors, age, and mass hyperspace. This 
work specifically is based on a My, {U - B)o, 
{B - V)o, age, and mass hyperspace due to the 
availability of integrated j7By colors and magni- 
tudes from [Hunter etani2003l . but the code can 
provide results based on more colors. The most 



^As part of the MASSCLEAN package, MASS- 
CLEANflgg is open source and freely avail- 
able under GNU General Public License at: 

^http : //www . physics . uc . edu/~popescu/massclean/ 

^X^ minimization provides reliable results in the case of Gaus- 
sian distribution, so it did not apply to low and medium mass 
clusters. 
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Fig. 7. — The result of 50 million Monte Carlo simulations, color coded by mass as described in Figureffl Here we see 
the full range of color-color magnitudes predicted as a function of mass. Clusters more massive than 50,000 M© (in the 
pink) typically adhere closely to the predicted values for the single solution Padova-based SSP predictions, computed 
in the infinite mass limit. However, clusters with masses below 10,000 M© (the dark green just beyond the yellow) can 
show extreme color ranges. 
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probable values for age and mass are reported, 
along with a list of other possible matches. The 
number of possible matches could be selected 
by the user or could be automatically computed 
based on the photometric errors. The results 
could be used to generate the confidence levels 
for age and mass (as presented in Tables 1 and 2 
and Figure [9j, or to generate the age-mass prob- 
ability distributions (as presented in FigurefTOl). 

By including additional bands (e.g. {V - K)o) 
the accuracy in the age and mass determination 
would increase, and could minimize the degen- 
eracies. In order to do it, a much larger version 
of the database is required, since the range of 
stochastic fluctuations is mu ch larger in {V - K)o 
l|Popescu & Hanson|[20T0ah . 



3.1. Age and Mass Determination for LMC 
Clusters 

As a first demonstration of the MASSCLEAN<3g^ 
code, we compare our values of ag es for 7 clus- 
ters w ith both spectroscopic ages dSantosetaJ] 
1200 6') and photometric ages (JHunter et al.ll2003h . 
Our results are computed using the previously de- 
reddened values_ofMi^^ - B)o, and {B - V)o 
from the [Hunter at al]l2003l catalog. The data is 
presented in Tab le [H The iritegrat ed photome- 
try and age from iHunter et al.l[2003l is presented 
in the columns 2-5, and the spectroscopic ages 



from iSantos et al. 2006 is presented in the col- 
umn 6. The MASSCLEANag^ results are pre- 
sented in the columns 7-9. The same data is 
also presented in Figure[9](a), (b). The integrated 
{U - B)o and {B - V)o colors are presented in 
Figure 8 as pink squares. 

The ages computed using MASSCLEANag^ 
are in goo d agreemerit with the spectroscopic 
ages from Santos et al.l |2006l . One object in Fig- 
ure[9]still lies signif icantly away from th e predicted 
age, as derived by lSantos et al .1120061 . even using 
MASSCLEANag^. We investigated a few of the 
sources in this diagram by looking closely at the 
minimization solutions found within the MASS- 
CLEAHage hyperspace. Absolutely every mass 
and age range available was specifically calcu- 
lated against the observed UBV colors. This 
output was used to create the likelihood plots 
shown in Figure [TOl These banana-shaped dia- 
grams demonstrate the very strong degeneracy 



in UBV colors between age and mass. Here 
we see that for NGC 1894, which represents the 
most significantly desperate point found in Fig- 
ure [9] (b), our best so lution for the age d oes not 
match that derived by lSantos et al .1120061 via inte- 
grated spectroscop y, but agree s with the age from 



Hunter et al. 200: 



opy, out agree s with the age Trom 
a If the older lSantos et al.l[2QQ6l 



age is correct, then this cluster is fairly massive 
(> 15,000 Mq based on Figure[iOj. 



3.2. Furtlier Tests of the MASSCLEANag^ 
Program 



Regrettably, the ISantos et ani2006l stud y has 
only a small overlap with the iHunter et ai]l2003l 
catalog of clusters, a consistent, high quality UBV 
dataset with SSP derived ages. But we wish 
to provide a further test of the MASSCLEANag^ 
package. The H unter et al. 2003 study includes 
over 900 LMC clusters, but only a fraction of those 
clusters are good examples. By this we mean, 
their location in the color-color diagram allows 
them to be more clearly aged using SSP models. 

We selected a subset of 30 clusters, which 
cover a wide range of {U - B)o and {B - V)o col- 
ors, but are also located close to the infinite mass 
limit (see the yellow triangles in Figure[8ll. The se- 
lected clusters are given in Table 2. Due to their 
locations in the color-color diagram, the age de- 
termination based on the classical x^ minimiza- 



tion methods and infinite m ass limit is expected 
to be r elatively accurate (e . g.Eancon & Mouhcinel 
I20OOI : iPopescu & Hanson"2010aV As shown in 
Figure Qj] our MASSCLEANag^ r esults are in 



goo d agreement with the ages from IHunter et al. 
^003 based on SSP model fits. However, we 
can also derive the mass, and this is presented 
in column 8 of the Table 2. New values for age 
and mass for all the clusters with UBV colors in 
the Hunt er et al. 2003 catalog and derived using 
MASSCLEANa^^ will be presented in a subse- 
quent paper (Popescu et al. 201 0, in preparation). 



4. Discussion & Conclusions 

That there would be degeneracies in the UBV 
color-color diagrams affecting the derivation of 
age as shown in Figures [7] and [8l was previ- 
ously recognized by many (e.a. lBruzual & Charlol 
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Fig. 8. — Monte Carlo cluster simulations colored by age for log{age/yr) = 7.00, 7.50, 8.00, 8.50, 9.00. Here we show 
that young clusters of low mass can be found to reside in the region that the infinite mass l imit predict for old stellar 
clusters. For a comparison with the observational data, the clusters fro m lHunter et al.ll2003l catalog are displayed as 
black dots. The clusters with photometric ages from lSantos et al.ll2006l (also presented in Table 1) are highlighted with 
pink squares. All the clusters presented in Table 2 are highlighted with yellow triangles. 
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Integrated Photometry'' 


Spectroscopy^ 




MASSCLEAN 






Name 


Mv 


(U-B)o 


(B - V)o 


Age 


Age 


Age 


Age 




Mass 




(mag) 


(mag) 


(m.ag) 


(Myr) 


(Myr) 


(Myr) 


(log) 




(Mo) 


1 


2 


3 


4 


5 


6 


7 


8 




9 


NGC1804 


-6.493 


-0.523 


-0.103 


19.1 ±7.21 


60 ±20 


29.5111^:^^3 


7 47+0.18 

'•^'-0.22 


2, 


000+2'°°° 

UUU_^ 3QQ 


NGC1839 


-7.087 


-0.347 


0.002 


16.2 ±7.16 


90 ±30 


72.44l-:^^3 


7 o^+O.lO 


8, 


000+2'°°° 

UUU_4 QQQ 


SL237 


-6.951 


-0.494 


0.156 


15.2 ±7.02 


40 ±20 


38.901^3:?^ 


7 P-q+O.lO 


5, 


P,PP,+i,ooo 

UUU_3 QQQ 


NGC1870 


-7.493 


-0.423 


0.082 


15.8 ±7.54 


60 ±30 


43.65tf,f 


7^4+0.19 


8, 


nnn+^'°°° 

^^^-2,000 


NGC1894 


-7.956 


-0.550 


0.192 


15.0 ±8.02 


100 ± 30 


20.89+1^4^9 


7 00+0.23 
''•^^_0.03 


7, 


P,^P,+5,000 


NGC1913 


-7.547 


-0.493 


0.257 


14.6 ± 7.72 


40 ±20 


33.11+^\f 


7 ^-9+0. 13 


7, 


P,^P,+2,000 
^UU_2,000 


NGC1943 


-7.255 


-0.359 


0.137 


15.5 ±7.31 


140 ± 60 


79AStl,% 


7 QO+0-04 

^•yu_Q 22 


10, 


QOO+i'°°° 



^Hunter et al. 2003 
^Santos et al. 2006 



20031 : ICervino & Valls-Gabaudl l2009l: iBuzzoni 
1 9891: IChiosilh 989l: lBruzualll2002l: lBruzualll201 Ol: 



The MASSCLEANc(9/(9r^ database is suffi- 
ciently populated now so it can be used to work 
back and derive mass and age for moderate and 
lower mass clusters in the LMC with a known 



Cervino & Luridiana 2004; Cervino & Luridiana 
2006 : Faaiolini et al. 20 07: Lan con & Mouhcind 
2000l ; lLancon & Fouesne au 2009: lFouesneau & LanconI (correctable) extinction. To make this possible 

2009'; Fouesneau & Lancon 201^;'Gonzalez et al.| with just three input bands, the dataset uses the 

2004 : Gonzalez-Lopezlira et al. 2005: Gonzalez-LopeziiraFMcM va 2008 models (Marigo et al. 2008) with a 



201 d: iMaiz Apellan zl l2009l : iPandev et ani201Ql : 
PopescLi& Hansonl l2010aV However, the scat- 
ter experienced in the integrated magnitude of a 
cluster also means using luminosity as a proxy for 
mass is highly dangerous for clusters of moder- 
ate and low mass. This is clearly demonstrated in 
our simulations, particularly in Figure HI Present 
day SSP models are not designed to work with 
mode r ate and low mass cluste rs (e.g. Bruzuali 
20021: lBruzualll2010l: iLancon & Mouhci nell2000l: 



Popesc u & HansonI l 2009l : [Popescu & HansonI 



2010a ). Yet, this is in fact the way in which mass 
functions for stellar clusters in other galaxies have 
been derived. For some studies, it might be ar- 
gued the masses are high enough that the errors 
are not extreme JLarsen & RichtleH l2000l ). but 
this method is being extended to a lower mass 
regime such as in the LMC (e.g. Billett t et a l.l 
20021 : lHunteretani2003l ). simply because there 
was no other way to derive this critical informa- 
tion. 



metallicity of Z = 0.008 ([Fe/H]=-0.6). While the 
simulations created for the database extend to 
ages up to log{age/yr) = 9.5, we have avoided 
applying the current inference package MASS- 
GLEANag^ on clusters in the LMC for which pre- 
vious indicators suggested ages greater than 
log{age/yr) > 9.0. This is because of our con- 
cern in the evolutionary models at such age and 
the inevitable drop in metallicity for the very old- 
est LMC clusters. It seems simplistic to as- 
sume the same metallicity can be applied to 
old clusters as to recently formed clusters, how- 
ever, the metallicity of the LMC does appears 
to have been reasonably slow changing over a 
fairly long period of time, from relati vely recent 



times until about loQiage/yr) = 9.5 JBica et c 
Eiii; [ Harris & ZaritsQ Iooi; | Piatti et al. |^0^. 



These studies also show that for some young 
clusters, log{age/yr) < 8, the adopted metallic- 
ity for the current MASSGLEANc(9/(9r^ database 
may be too low by a few dex. Whether the 
MASSCLEANco/6>r^ database might be extended 
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Fig. 9. — In (a) above, the log{age/yr) as derived 
by Hunter et al. (2003) using UBV photometry alone 
is show n versus the log{ag e/yr) as independently de- 
rived bv Santos et al.l2006l based on integrated spectra 
of those same clusters. In (b) below, we show the same 
set of stellar clusters with ages as derived by MASS- 
C LEANggg and using the same input UBV photometry 
as lHunteret"ani2003l used. 
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Fig. 10.— Output from MASSCLEANag^: finding 
probabilistic regions in the mass - age plane, based on 
observed UBV magnitudes for these clusters. These 
diagrams show the very strong correlation between 
cluster mass and age in integrated UBV colors. Our 
best sol ution for NGC 1894 does not match the age 
given bv lSantos et"ai1Eooa 
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Fig. 11. — We have further checked the consistency of 
the ages de rived using MASSC LEANag^es by selecting 
a subset of iHunter et al.ll2003l clusters with previously 
well constrained ages. These data are also shown in 
Table 2. 



to include LMC clusters over a broader range 
of metallicities, from [Fe/H] = -1 to and for 
log{age/yr) > 9.5, and if this might prove signifi- 
cant in deriving accurate age and mass estimates 
remains to be explored. If such an additional pa- 
rameter would be added to the database, more 
input photometry will be needed for the stellar 
cluster in question, to keep the search from be- 
coming hopelessly degenerate over all the prop- 
erties. 
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